arXiv:1503.07784v2 [cond-mat.mes-hall] 6 Oct 2015 


Odd-Prequency Triplet Superconductivity at the Helical Edge of a Toplogical Insulator 


Frangois Crepin, Pablo Burset, and Bjorn Trauzettel 
Institute for Theoretical Physics and Astrophysics, 

University of Wurzburg, 97074 Wurzburg, Germany 
(Dated: October 7, 2015) 

Non-local pairing processes at the edge of a two-dimensional topological insulator in proximity 
to an s-wave superconductor are usually suppressed by helicity. However, additional proximity of a 
ferromagnetic insulator can substantially influence the helical constraint and therefore open a new 
conduction channel by allowing for crossed Andreev reflection (CAR) processes. We show a one- 
to-one correspondence between CAR and the emergence of odd-frequency triplet superconductivity. 
Hence, non-local transport experiments that identify CAR in helical liquids yield smoking-gun ev¬ 
idence for unconventional superconductivity. Interestingly, we identify a setup - composed of a 
superconductor flanked by two ferromagnetic insulators - that allows us to favor CAR over electron 
cotunneling which is known to be a difficult but essential task to be able to measure CAR. 
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Introduction .— The influence of strong spin-orbit cou¬ 
pling and the constraint of time reversal symmetry are 
responsible for the appearance of helical electronic chan¬ 
nels at the edge of two-dimensional (2D) topological in¬ 
sulators [1-3]. Ample evidence for the experimental de¬ 
tection of these edge states has been seen in transport 
[4-6] and scanning SQUID experiments [7]. Proximity 
of an ordinary s-wave superconductor can induce pair¬ 
ing at the helical edge [8, 9]. Interestingly, the interplay 
of helicity and superconducting order gives rise to un¬ 
conventional proximity-induced superconductivity (SC) 
in these systems. However, it is fair to say that it is 
very difficult to unambiguously probe the emergence of 
unconventional SC because - often times - conventional 
and unconventional signatures of SC look alike. Several 
recent experiments have demonstrated (as a first step to¬ 
wards the detection of unconventional SC) that helical 
liquids as boundary states of quantum spin Hall systems 
can indeed be brought in proximity to s-wave supercon¬ 
ductors [10] and serve, for instance, as conducting chan¬ 
nels of a Josephson junction [11, 12]. 

Importantly, helicity guarantees perfect local Andreev 
reflection [13] - the conversion of an electron into a hole 
with opposite spin - at the interface between the normal 
and the proximity-induced superconducting region called 
NS junction. Evidently, perfect local Andreev reflection 
can give rise to sub-gap Andreev states, for instance, in 
Josephson junction setups. Among these sub-gap states, 
the one with zero excitation energy is its own charge- 
conjugate state. It has been coined Majorana (bound) 
state and has attracted a lot of attention because of po¬ 
tential applications of this anyonic state for topological 
quantum computing [14]. If a ferromagnet (FM) is placed 
in vicinity to the NS interface the Majorana (bound) 
state can be localized in the region between the FM and 
the SC [9]. This localization allows us to make a formal 
connection to the physics of a finite-size ID p-wave topo¬ 
logical superconductor [15], and its potential realizations 


in spin-orbit nanowires [16, 17]. 
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FIG. 1. (Color online) Examples of hybrid junctions in a 
helical liquid. A possible scattering process is represented, 
with a spin up right moving electron (solid red line), a spin 
down left moving electron (solid blue line), a spin down right 
moving hole (red dashed line) and a spin up left moving hole 
(blue dashed line). 


In this Letter, we give for the first time a complete 
characterization of the pairing symmetries in a helical 
liquid, in proximity to both FM insulators and SC. The 
typical setups we have in mind are sketched in Fig. 1. 
Interestingly, we discover that the presence of both FM 
and SC yields an unconventional contribution to the pair¬ 
ing amplitude that is odd (O) in frequency, triplet (T) 
in spin space, and even (E) in orbital space (in short: 
OTE). Now, the natural question arises whether it is 
possible to unambiguously detect this exotic pairing in 
our setup. Remarkably, the answer is yes. Indeed, the 
FM also allows for equal-spin crossed Andreev reflec¬ 
tions, that would otherwise be totally suppressed by he¬ 
licity. Conversely, in the absence of the FM, equal spin 
OTE pairing is forbidden by time-reversal symmetry and 
CAR contributions completely vanish. Our proposal, 
based on well established measurements of the differential 
conductance [18, 19], turns out to provide the simplest 
constituents necessary to directly probe odd-frequency 
triplet pairing. Hence, CAR is a crystal-clear transport 
signature of unconventional superconductivity in helical 
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liquids. Our system substantially differs from the origi¬ 
nal literature about non-helical, diffusive, FSF junctions, 
where long-range proximity effect was considered as evi¬ 
dence for odd-frequency triplet superconductivity [ 20 ]. 

Previously, the shot noise related to CAR processes 
has been proposed as a way to detect Majorana bound 
states [21, 22]. Moreover, in a nanowire setup, local An¬ 
dreev reflection measurements have been suggested to 
probe spin-triplet superconducting correlations of Majo¬ 
rana fermions [23]. However these works both address 
a fundamentally different scientific question than ours. 
In our Letter, we describe the novel idea of relating 
CAR processes directly to symmetry properties of the 
proximity-induced pairing amplitude. This idea is heav¬ 
ily based on the peculiar transport properties of helical 
edge states. 

Model and analysis of pairing symmetries .— Our start¬ 
ing point is the Hamiltonian (we set h = 1) 


incoming from the right), (pui.i (hole-like quasi-particle 
incoming from the right). Each incoming particle can ei¬ 
ther be Andreev reflected, with probability amplitude a^, 
or normal reflected, with probability amplitude bi. We 
give the asymptotic forms of the scattering states in the 
supplementary material [26] . We proceed now to analyze 
the symmetry content of the pairing in hybrid junctions 
based on helical liquids. The Green’s function is a 4 x 4 
matrix of the general form 


G^{x,x' .uj) 


fG^[x,x',uj) Gfj^{x,x',uj)\ 
\Ghei^,x',uj) Gj^f^{x,x',uj)J 


( 4 ) 


The pairing amplitude, given by the off-diagonal elements 
of Eq. (4) must be completely antisymmetric under the 
exchange of positions, times and spin components of the 
electrons making up the Cooper pair. The spin symmetry 
can be made explicit by writing 


1 

H = - dx {x)'HBdG{x)'it{x) , (1) 

-J — OO 

where we have introduced the spinor 4 /( 2 ;) = 

['tf^{x),'tjji{x),if]^{x),—if^{x)]'^ and the Bogoliubov-de 
Gennes (BdG) Hamiltonian HBdGix) = {—ivpdj^az — 
P.)tz + m(a;) • cr + A{x)tx. The Pauli matrices cr = 
(fx, CTy, CTz) and T = {Tx,Ty,Tz) act in spin and particle- 
hole space, respectively. The phase of the pair¬ 
ing potential plays no role in our calculation and we 
have omitted it for the sake of clarity. In prac¬ 
tice we take A(x) = 6{x)9{d — x)Ao and m(a;) = 
mi(x) -I- m 2 (a;) with mi(x) = 0(xi + S — x)9(x — xi + 
(5)(toi cos Ai, TOi sin Ai, Wzp) and m 2 (a;) = 9{x — X 2 + 
(5)0(x 2-I-5 — a;)(m 2 cos A 2 , m 2 sin A 2 , 7712 , 2 )- In Fig. 1(a), 
we illustrate the simpler case where m 2 = mz ,2 = 0, 
whereas Fig. 1(b) corresponds to the most general case 
we consider here. We are interested in the pairing ampli¬ 
tude, that is, the anomalous (electron-hole) part of the 
retarded Green’s function 

G^{r,r') = -i9{t-t'){{'i/{r),'it\r')}) , (2) 

where r = {x,t) and r' = (x',t'). Its Fourier transform 
G^{x,x',uj) = Jdt satisfies the 

equation of motion [oj — HBdG(a^)] G^{x, x', w) =^(a:—a;') 
and is discontinuous in a; = a;', such that 

lim [G''^(a;' + s,x',uj) — G^{x' — £,x', w)] = , ( 3 ) 

£->■0 IVp 

For pedagogical reasons, we first investigate the setup 
of Fig. 1(a) in the limit d —>■ -too. We construct the 
retarded Green’s function from scattering states satisfy¬ 
ing outgoing wave boundary conditions [24, 25]. These 
states define an eigenbasis of the BdG Hamiltonian. For 
a given energy w, there are four independent scattering 
states (electron incoming from the left), (hole 
incoming from the left), (electron-like quasi-particle 


Geh{x,x',uj) = f^{x,x',uj)ao + f^{,x,x',uj)a^ , (5) 

where is the singlet (S) component of the pairing am¬ 
plitude while f [^2 3 make for the three components of the 
triplet (T). The antisymmetry of the pairing amplitude 
translates as 


fd{x,x,w) = f^{x,x,-uj) , ( 6 ) 

fl,2,di.XiX',L0) = -f^2,3i^',x,-uj) , (7) 

where /o^i, 2,3 is built from the advanced Green’s function 
G^[x,x',uj) = G^{x',x,ujy. The rule of thumb is gen¬ 
erally that the orbital part of the singlet (resp. triplet) 
is even (resp. odd). However, this is not the most gen¬ 
eral case, as the behavior under the exchange of times 
also contribute to the antisymmetry of the Gooper pair 
Green’s function [27] . In principle, each of the f^oi 2 3 
can be even (E) or odd (O) under the exchange of x 
and x', that is, f^{x,x'.u) = ±fj^(x',x,oj), or under 
w —>■ —OJ, that is, f^{x,x',oj) = ±f^(x,x',—uj). The 
pairing is thus classified according to these three sym¬ 
metries and, adopting the notation [frequency, spin, or¬ 
bital] , there are only 4 allowed combinations that comply 
with the total antisymmetry of the pair amplitude: (ESE, 
OSO, ETO, OTE) [28]. Actually, in the case at hand of 
an NS junction, translational invariance is broken and 
the 2 3 are necessarily a mixture of functions with 

even and odd orbital symmetries. We find that the com¬ 
ponents of the pairing amplitude are of the form [29] 


fo^^ix,x',uj) = 
ff^'^{x,x',oj) = 
f±^{x,x' ,oj) = 

where /+ = /i - 
ply the tt (resp. 


/o^/„terface(a;. a:', Oj) + foityi,{x, x', Uj), 

/3^/„terface(a;, x', Oj) ^ x', Oj), 

/Snterface(a;,a;',a;), (8) 

- i /2 (resp. /_ = /i -I- i/ 2 ) is sim- 
//) pairing amplitude. The bulk 
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contributions are all proportional to e ^ I with 

k{uj) = — oj'^ Ivp, while the interface contributions 

decay to zero as ) [26]. Although the pairing 

term in the Hamiltonian (1) is of the singlet form and re¬ 
flects the symmetry of the nearby superconductor [30] , a 
triplet component naturally emerges in the bulk, because 
of helicity. Deep inside the proximity induced supercon¬ 
ducting region, we find a coexistence of a singlet and a 
triplet component. We have verified that they belong 
to the conventional symmetry classes, namely ESE and 
ETO, respectively [31]. At the interface however, all of 
the symmetry classes matter. The singlet has both ESE 
and OSO components while the triplet /a has both ETO 
and OTE components [32] . However, the two other com¬ 
ponents of the triplet, /+ and /_, corresponding to ft 
and pairing respectively, exist only in the exotic OTE 
class. Explicitly, we find 


/^.i„terface(a;,a;',w) 

/—,interface{^7 ^ 7^) 
/-t-,interface(^7 ^ J 

/Ainterface(a;,a;',w) 


= —Aq 
= —Ao 
= —Ao 
= —Ao 


k(ui))(x+x' ) 

2vf\/^o ~ 

k{ui)){x+x' ) 

2 u_Fy/Ag — 

k{ui))(x+x' ) 

2vf\/^o ~ 

k(ui)){x+x' ) 

2vf\/^o ~ 


«&3(w) , (9) 

-&4(w) , ( 10 ) 

M^r 7 ( 11 ) 
M^r 7 ( 12 ) 


with 63 / 4 ( 0 ;) the normal reflection amplitudes for 
quasi-particles incoming from the right [26]. It is clear 
that ,u}) = f^^^{x',x,oj). Furthermore, 

the fundamental particle-hole symmetry of the BdG 
Hamiltonian, Uc'H’^^qUq = —'H'bag, with Uc = cryTy, 
yields 63 ( 0 ;) = — 64 (—o;)*. It hence follows that 
f^{x,x',oj) = —f^{x,x',—u}). In the absence of the 
ferromagnetic insulator, normal reflection processes are 
forbidden by helicity and 63/4 —?> 0. Therefore, proximity 
of a ferromagnet induces odd-frequency triplet pairing in 
the tt and H channels, at the edge of a 2D topological 
insulator. This constitutes one of the main findings of 
this Letter. The symmetry classification is summarized 
in Table I. 


For simplicity, we now consider the case of a ferro¬ 
magnetic impurity, meaning that we reduce the width 
of the ferromagnetic domain d —>■ 0 while keeping the 
dimensionless quantity {miS)/vF = rno finite. We take 
rriz^i = 0, with no consequence on the physics [33]. In 
that case, we expect the formation of Andreev bound 
states between the impurity and the superconducting re¬ 
gion, the number of which is essentially given by the ra¬ 
tio of the distance to the superconductor jxij and the 
coherence length ^ = vp/^o- Their main characteristic 
is to allow for perfect local Andreev reflection [34]. In¬ 
terestingly, it is possible to relate the physics of these 
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Interface 

Bnlk 

fo 

tt-tt 

ESE+OSO 

ESE 

h 

tt + tt 

ETO-bOTE 

ETO 

f± 

tt,tt 

OTE 
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TABLE I. Symmetry classification of the pairing in a NFS 
junction at the edge of 2D topological insulator. E and O 
indicate even and odd, while S and T stand for singlet and 
triplet. The first letter gives the symmetry with respect to 
frequency while the last gives the orbital symmetry of the 
pairing. 



FIG. 2. (Color online) Pairing amplitudes fo{x,x,u}) (blue 
dotted line), f 3 {x,x,u}) (red dashed line) and f+{x,x,uj) 
(black solid line), at the NS interface a; = 0 [see Fig. 1(a)]. The 
impurity is located at xi = —3^ and its strength is mo = 0.5. 
The positions of the peaks correspond to the localized An¬ 
dreev states in the junction and signal an enhancement of 
triplet pairing. 


localized states to the induced triplet superconductivity. 
In Fig. 2, we plot the moduli of the local pairing ampli¬ 
tudes fo{x, X, w), fz[x, X, w), and /+(x, x, w) as a function 
of the excitation energy w, for x = 0 (see Fig. 1(a)). Note 
that /+(x, x', w) = /_ (x, x',uj) in the impurity limit. En¬ 
forcing X = x' effectively gets rid of the odd orbital com¬ 
ponents. Therefore we probe only ESE and OTE compo¬ 
nents of the pairing. We have taken jxij = 3^ and a quite 
large impurity strength mo = 0.5, in order to have three 
well-defined bound states below the gap. The ft pairing 
amplitude is especially strong close to the energies of the 
bound states. It is actually maximal at the NS interface 
X = 0 , and decays exponentially into the bulk as can 
be seen in Eq. (9). This confirms the proximity induced 
character of the odd-frequency triplet pairing in a helical 
NFS junction. It should be noted however that no local 
observable can directly probe the ft or || pairing. We 
have indeed checked that the ft and local pairing am¬ 
plitudes vanish at the NF interface and, by continuity, in 
all of the normal lead: f^{x,x,uj) = 0 for x < Xi. Fur¬ 
thermore, we obtain f^{x, x, oj) = f^{x, x, uj) for x < xi, 
that is equal — it and ft + tt correlations in the nor- 
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Vi/Ao 

FIG. 3. (Color online) Non-local conductance G 21 = 

—dhIdVi = Tcar — Tec, at zero temperature, in the setup 
of Figs. 1(a) (blue dashed line) and (b) (red solid line). A 
positive signal indicates an excess of crossed Andreev reflec¬ 
tion over electron tunneling. We have taken d = 2^, 5 = 0.5^, 
mi = I. 5 A 0 , mzp = 0, a:i = —2^, X 2 = d -f and a non¬ 
zero chemical potential /r = O.SAq. In the first setup (blue 
dashed line) there is no ferromagnetic impurity in the right 
lead, m 2 = mz ,2 = 0, while in the second setup (red solid 
line) m 2 = mi and mz ,2 = mz,i = 0. The black dotted line 
shows the conductance in the absence of any ferromagnetic 
impurity. 


mal lead, again a direct consequence of helicity. Thus, 
the local conductance Gii(Vi) = dli/dVi, with /i the 
current injected in the left (normal) lead and Vi the bias 
between this normal lead and the superconductor, does 
neither probe ft nor H correlations. In order to detect 
these correlations, one needs to consider the non-local 
conductance. 

Crossed Andreev reflection and transport signatures .— 
We now consider the setup of Fig. 1(a) with d finite. 
We have seen that helicity forbids to probe the ft and 
44 . pair correlations locally. Indeed, a right moving elec¬ 
tron incoming from the left lead cannot be reflected as 
left moving hole with the same spin. However, the right 
moving electron can be converted into a right moving hole 
with the same spin, in the right lead. The latter process, 
called crossed Andreev reflection (CAR), corresponds to 
the injection of a Cooper pair in the ft channel [35, 36], 
and probes directly the exotic pairing induced by the 
proximity of the FM insulator. We do not expect that 
charging effects are important in our setup. The reason 
is that for a size of the SC region of order 4, tbe cor¬ 
responding single particle level spacing should be much 
larger than the corresponding charging energy. Thus, we 
have solved the problem based on scattering theory and 
computed the CAR and electron cotunneling (EC) prob¬ 
abilities. As can be inferred from Eq. (9), both EC and 
CAR probabilities are suppressed exponentially with the 
length d of the superconducting region. However CAR 
is enhanced for excitation energies close to that of the 
Andreev bound-states, consistent with the enhancement 


of triplet pairing shown in Eig. 2. Unfortunately, EC is 
always bigger than CAR in that case which is a rather 
generic situation to arise. The non-local conductance 
G 21 = — 9 / 2 / 9 V 1 = Tcar — Tec directly probes the dif¬ 
ference between CAR and EC. It can be seen in Eig. 3 
(blue dashed line) to be always negative for this setup. 
Contrary to other proposals, we cannot filter spin in or¬ 
der to favor CAR [37]. Another option is to break the 
symmetry between particle and hole excitations [38], yet 
we find that the chemical potential has virtually no effect 
on the transmission probabilities in the case of a single 
impurity. 

However, there is a very convenient way out of this 
dilemma: Adding a second EM to the right lead (setup 
of Eig. 1(b)) allows to significantly enhance CAR over 
EC. Andreev bound-states will now also arise in the right 
lead, their number being given by (x 2 — d)/^. If the 
length of the superconductor is not too big compared 
to 4, the two sets of bound states (left and right of the 
SC) hybridize, and become very sensitive to the chem¬ 
ical potential. We found that CAR and electron tun¬ 
neling can then be made quite different in systems with 
a pronounced left-right asymmetry. As an example, we 
have taken Xi = —2^ and X 2 — d = ^. The zero bias 
peak in the conductance splits under the effect of the 
hybridization, as can be seen in Eig. 3. However, the 
second bound-state, that is localized in the left lead, 
has no counterpart to hybridize with. This results in 
a single peak in the non-local conductance, whose sign 
can be tuned at will with the bias. A positive signal of 
G 21 = — 9 / 2 / 5 V 1 = Tcar — Tec is an unambiguous sig¬ 
nature of odd-frequency triplet pairing at the edge of 2D 
topological superconductors. 

Conclusions. — Helicity at the edge of a quantum spin 
Hall system implies perfect spin-to-charge conversion: all 
right movers have spin up and all left movers spin down. 
In proximity to an s-wave superconductor, helical edge 
states exhibit unconventional pairing amplitudes. We 
have carefully characterized the symmetry of the super¬ 
conducting pairing at the helical edge in the additional 
presence of ferromagnetic insulators. Interestingly, we 
have discovered a unique relation between the existence 
of odd-frequency triplet superconductivity in these sys¬ 
tems and crossed Andreev reflection processes due to he¬ 
licity. Remarkably, we have even identified a possibility 
to favor crossed Andreev reflection over electron cotun- 
neling by a proper tuning of resonance conditions. 
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Supplementary material 

In this supplementary material we give detailed expressions for the (retarded) Green’s function. We start by 
discussing the structure of the Green’s function in particle-hole and spin space. We then give explicit formulas for 
the junction of Fig. 1(a) in the limit d —>■ oo, and give additional information about the pairing symmetries, as well 
as analytic expressions for the scattering matrix in the case of a ferromagnetic impurity. 


A. General structure of the Green’s function 

It should be noted that because of the particle-hole symmetry of the BdG Hamiltonian, not all sixteen elements of 
the Green’s function are independent. Indeed we have, 

UcH*b^gUI = -HBdG , 


(13) 
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and 


UcG^{x,x',uj)*ul; = -G^{x,x',-uj) , 


with Uc = CTyTj,. Defining, 


we find that 


G^{x,x\uj) 


(g^{x,x',uj) G^^^{x,x',uj)\ 
\GLix,x',uj) Gf^f^{x,x',uj)J 


Gee(x,x',-uj) = -ayGj^f^{x,x',u})*ay , 

Gehix,x',-uj) = ayG^^{x,x',uj)*ay . 


(14) 


(15) 


(16) 

(17) 


Therefore, the knowledge of e.g. G^^{x,x',oj) and G^f^{x, x',u}) is enough to characterize the full Green’s function, 
and all the associated physical processes. We construct the Green’s function from the scattering states with outgoing 
boundary conditions. For the sake of completeness we give their asymptotic forms, in the leftmost normal region 
{x < xi — S) and in the superconducting region (x > 0). For a given energy w, there are four independent scattering 


states 4’2,u:,4>3,uj, 4>4:,ui- normal region, 

, (18) 

Mx) = 4^^+ 024-^6-*'='" , (19) 

, ( 20 ) 

AM) = , ( 21 ) 


where we have introduced kf,ffi = (p ± u})/'>^f, as well as (jjA = (1,0,0,0)^, 4 ^ = (0,1,0,0)^, = (0,0, 0,1)^, 

and (j)\^ ^ = ( 0 , 0 , 1 , 0 )^. In the superconducting region we have 

AM) = cixAeA- + , ( 22 ) 

AM) = d2xAe^''"- + €2x44-*'="" , (23) 

AM) = xAe-A- + 03x4^6-*'="" + 63X44*'="" (24) 

AM) = xi~^e-A- + 64x44-*""" + 04X44*'="" (25) 


where now = ik ^ — A'^)/vf, for excitation energies above the gap. 


X. 


(+) 


= (e^(-), 0 , 1 , 0 )^, X 


(-) _ 


= ( 0 ,e^(“), 0 , 1 )^ 


x 4 ^ 


We also have introduced 
= ( 0 ,e-^(“), 0 ,l)^. 


xi ^ 


MMy ’ ’ ’ ' ’ ’ '^'* MM)' 

, ^ ( 6 "^^*^^ 0 , 1 , 0 )^, with = 1 + as well as AlA = arcosh(w/Ao). At a later stage, these 

expressions will be analytically continued to energies below the gap. The coefficients (a^, 6 ^, Ci, (ii)i=i, 2 , 3,4 in the scat¬ 
tering states depend on the excitation energy oj and are the elements of the scattering matrix. They are obtained by 
matching the wave-functions at the boundaries between regions. We demand that G^(a;,x*,a;) satisfies outgoing wave 
boundary conditions and construct it in the following way: 


pRirr p ) AAAsixA -k AAAix')'^ if a: < x' 

I ^i(a;)Ai(a::')'^-I-(()2(a:)A2(a;')'^ if x > x' 

The 16 unknowns in the vectors Ai, A 2 , A 3 , A 4 are solutions of Eq. (3), and are related, in the end, to the elements of 
the scattering matrix. We find G^(x,x',w) and Gf)j(x, x', w) in the various regions of the setup pictured in Fig. 1(a). 
Given our choice of the spinor 'I'(x) = [A{x),tpi{x)Al{x), ~V’j(^c)]'^, the spin components of electron and holes are 
distributed in the following way: 


\[G'fe]fT [GAnj 


(26) 











for the regular part, and 
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[{Gl]u {GlUtJ ■ 

for the anomalous part. 


(27) 


B. Green’s function in the left lead 


In the normal lead (x < Xi — S), 


lG^J^^(oj,x,x') = - x') , (28) 

[G« ];;(a;, X, x') = e-^^^G-^'^e{x' - x) , (29) 

[G^,]n{uJ,x,x')=0, (30) 

[G«U^(.;,x,x') = e-'=»("+“')6i(cu), (31) 

while 

= 0 , (32) 

[G^^Mio, X, x') = e*('=''"'-'='")a2(a;) , (33) 

[G^hhti^’ x,x ) = 0 , (34) 

[Gfh]u(w,2;,a;') = 0 , (35) 


where ± u})/vf, &i(w) is the normal reflection amplitude for electrons and 02 (w) is the Andreev reflection 

amplitude for holes. The Heaviside functions in the diagonal elements of G^^{x,x',uj) arise from the discontinuity 
of the Green’s function at a; = a;' given in Eq. (3), a direct consequence of dealing with Dirac fermions. Note that 
some matrix elements in G^^{x,x',uj) and G^(a;, a;', w) vanish because of helicity, and can never be non-zero, even 
in the presence of the ferromagnetic impurity. At this level we can already discuss the spin symmetry of the pairing 
amplitude by writing 


Gfhix,x,u}) = f^{x,x,u})ao + f’^ix,x,uj)(7i . (36) 

It follows directly from the structure of Gff^{x, x',uj) that fi{x,x',uj) = f 2 {x,x',u}) = 0 and fo{x,x',uj) = —/^{x, x', oj) 
for all values of x and x' in the left lead. From this we deduce that ft and H cannot be probed locally, in a transport 
experiment, and that local Andreev reflection in a helical liquid induces equal singlet (fi — it) and triplet (ft + it) 
correlations. 


C. Green’s function in the superconductor 


In the region covered by the superconductor {x > 0) the physics is much richer, and we find, for excitation energies 
below the gap 


[Gee\uixi,X,x') 

[Gee]ni.^^^^^') 

[Gee]\.^{^^^X,x') 


2 ^ ^iti{x-x') 


ivp 

e‘ 

iri{Lo 

') _ ^-i'nG) 

I 


e~ 


ivp 

e‘ 

i7)(uj 

•) — g-jr/G) 

I 


e 

i^{x-\-x') 

ivp 

e‘ 

i7f(uj 

•) — 

I 


e~ 

-ifi{x-\-x') 

ivp 

e‘ 

irf(cj 

') _ ^-iriG) 


e-<‘^)\x-x'\ + e{x' - a;)e-*''(“) 

e-<‘^)\x-x'\ (^0(a: - + e{x' - 

h3{uj)e-<^'^G+=^''> , 




(37) 

(38) 


(39) 


(40) 












for the regular part, and 




[Geh\i^i^,X,x') 


[Geh]n{^^X,x) 


[Geh\n{^^X,x') 


1 




-x') 

IV F 

e‘ 

ir]{u 

) - 

g-ir)(in) 

1 


e~ 

i/i(: 

x — x') 

ivp 

e‘ 

ir]{u 

) - 

g-i7)(in) 

1 




ivp 

e‘ 

ir]{ui' 

) - 

g-i7)((n) 

1 


e~ 

ifi{xFx') 

ivp 

e‘ 


) - 



^-k{ui)\x-x'\ _|_ Qg(^)g-*'?(“)g-«(‘^)(®+2:') 


(41) 

(42) 

(43) 

(44) 


for the anomalous part, with k{oj) = a/A g — and 77 ( 0 ;) = arccos(a;/Ao). Each of the matrix elements is the sum 
of a bulk term (proportional to I) ^nd an interface term (proportional to g-«(‘^)(®+ 2 ^ )). in order to discuss 

the symmetries of the pairing amplitude, we also give the anomalous part of the advanced Green’s function 


[Gehlni^^^^^') 

[Gehhti^^X^x') 

[Geh]n{^,X,x') 

[Geh]u{^.X,x) 


2 ^iuix-x') 


ivp 

e‘ 

ir](uj) 

_ g-i? 7 (in) 

1 


e~^ 

,fl{x — x') 

ivp 

e‘ 

iri(uj) 

_ g-i 7 )(a;) 

1 


^iFix+x') 

ivp 

e‘ 

ir](uj) 

_ g-iJ)(a;) 

1 


e“^ 

Fix+x') 

ivp 

e‘ 

iri{uj) 



+ 04 ( 0 ;)* 

+ 03 ( 0 ;)* 

63 (w)*e"'=(“)(“+’"') . 




^-k(uj)(x+x') 


(45) 

(46) 

(47) 

(48) 


We expect the pairing amplitude (the Cooper pair wave-function) to be totally antisymmetric under exchange of 
positions, spins and times of the two electrons making up the pair. More precisely 

[Gf/,]n(w: X, x') = -[G^^]tt(-a;, x', x) , (49) 

[Geh\n(^^ X, x') = -[G±]^^{-uj, x', x) , (50) 

[Geh\ni^^ ^') = [Geh\t\{-^^ x', x) . (51) 


Note that the last equality comes without a minus sign because of the choice of basis we have made. The spinor 
defining the Hamiltonian of Eq. (1) indeed reads 4'(a;) = ['tfj^{x),ijji{x),tfj'^^{x),—tjj'^^{x)]'^. The minus sign appearing 
in the definition of the spinor implicitly takes care of the minus sign that should arise upon exchanging f and f 
in Eq. (51). From particle-hole symmetry, we have 04 ( 0 ;) = —a 3 (—w)*, 64 ( 0 ;) = and with the help of 

g*’ 7 (-‘^) = _g-*’)(<^) one can check the general validity of Eqs. (49) to (51). In order to isolate the spin symmetry we 
use the decomposition of Eq. (36). We have, on the retarded side. 


/o^buik(a;,a:',a;) 

/o,interface(^> 2 ; , w) 

/3^buik(2;,a:',a;) 
/3,interface{^5 ^ 5 ^) 
/-t-,interface (^5 ^ 5^) 
/—.interface(^5 ^ 5^) 


Aq 

2vf\/^q — 

Aq 

duFt/Ag - a;2 
»Ao 

2vf\I Aq — LtP- 

Aq _ 

4uF-\/Ag - a;2 
_ Aq _ 

2vf\/ Aq — uj^ 
_Aq_ 

2 vf\/ Aq — uj^ 


cos[n{x - x')] , 


g-K(in)(a:-|-a;') 


gr^r(x-x')^ 4 (^)g*r,(cn) 


sm[ii{x - x')] 


^-k(ui)(x+x') 


giAr(x-a.')^ 4 (^)g*r,(cn) 




^-^(u>){x+x’)^-^F{x+x')l,^(^^'^ ^ 






(52) 

(53) 

(54) 

(55) 

(56) 

(57) 
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where /+ = /i - i /2 (resp. /_ 
/o,bulk(^> X ,U!) = 
/o,interface(^> ^ 5^) 
ftbnlkix^x',Uj) = 

/3^i„terface(a;,a;',w) = 
/+,interface(a;,a;',w) = 
/-,interface(a;,a;',w) = 


= /i + i/2) is simply the ft (resp. //) pairing amplitude. On the advanced side, 
_ ^0 _ 

2vpy^ Ag — 

^0 


g-K(a;)|x-a:'| _ 3 ./)] 




K(cL>)(a:+ic') 




iAo 


2vf\/ Aq — 

_Aq_ 

4u_f^A^ - u ;2 
Ao 

2vf^ Aq — up- 
Ao 

2 Ui?-\/Ag — 


g-K(a;)|x-a:'| sin[^(x - a;')] , 




,i/i(a;-a;')Q*^g_,^gir;(ca) _ g-v(a:-a:') i,, 


a;(w)e- 


g-K(i,;)(a;+a;')giM(a:+a:')^*^^;) ^ 


( 58 ) 

( 59 ) 

( 60 ) 
( 61 ) 
( 62 ) 
( 63 ) 


One can check on these expressions some of the affirmations made in the main text. In particular, the bulk contribution 
in /o is of the ESE type = fi^huiki^' ’^hile , x', uj) = fo^h^ii,{x',x,-uj)), and the bulk 

contribution is of the ETO type (/3^buik(3;>2;',w) = -f^^^^^^^{x',x,UJ), while /3^buik(3^.= /3^buik(3;'.-w)). 
Eurthermore, /o,interface gives rise to both an ESE and an OSO contribution, while /s,interface gives rise to both an 
OTE and an ETO contribution. 


D. Case of a ferromagnetic impurity 

So far, all expressions are completely general and we have not specified the elements of the scattering matrix. 
However, simple analytic expressions can be obtained in the case of a ferromagnetic impurity. We let 5 —)■ 0 and keep 
miS = mo finite. We find, below the gap 


&i(w) 


ai(w) 


hioj) 


aoioj) 


1 sinh(4mo)(e-2i|a:il^gir,M g2r|a:i|a;g-rr,(a;))grAi 

2 i g2r|a;i|ccg-i7)(cc) sinh^(2mo) + cosh^(2mo) ' 

1 

g2i|a;i|ajg-r7)(cc) sinh^(2mo) + e“2*kicosh^(2mo) ’ 

1 sinh(4mo)(e®''^‘^) — 

2 i g 2 r|a;i|ccg-i 7 )(cc) sinh^(2mo) + g-2*ki|wgrr;(cc) cosh^(2mo) ’ 

g2*l^il“ sinh^(2mo) + g“ 2 *|a;i|<.c cosh^(2TOo) 

g2i|a;i|a;g-f»7(a;) sinh^(2mo) + g-2f kl |ccgf»7(a;) cOSh^(2TOo) ’ 


( 64 ) 

( 65 ) 

( 66 ) 

( 67 ) 


The other elements are given by particle-hole symmetry: 04(02) = —03(—w)*, 64(0;) = —bo(—uj)*, 02(0;) = —ai(—0;)*, 
62(02) = — 5 i(—02)*. The normal reflection coefficient vanishes for energies given by 


4 |xi|a; — 277(02) = TT -\- 2 mr, n G Z . 


( 68 ) 


Those are the energies of the bound-states that form between the ferromagnetic impurity and the superconductor. 
They induce perfect local Andreev reflection. 



